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Spin dynamics evolution of BaFe2(Asi_iPi)2 was probed as a function of P concentration via 
^^P NMR. Our NMR study reveals that two-dimensional antiferromagnetic (AF) fluctuations are 
notably enhanced with little change in static susceptibility on approaching the AF phase from the 
superconducting dome. Moreover, magnetically ordered temperature 6 deduced from the relaxation 
rate vanishes at optimal doping. These results provide clear-cut evidence for a quantum-critical 
point (QCP), suggesting that the AF fluctuations associated with the QCP play a central role in 
the high- Tc superconductivity. 



PACS numbers: 74.70.Xa, 74.40. Kb, 74.25. nj 



Conventional phonon-mcdiatcd superconductivity oc- 
curs in a normal metal that is well accounted for by Lan- 
dau's Fcrmi-liquid (FL) theory. However, the standard 
FL theory appears to break down above Tc in many "ex- 
otic" superconductors characterized by unconventional 
pairing rather than a conventional uniform-sign s-wave 
pairing function, such as in the heavy fermion materi- 
als and cuprates. The origin of the anomalous metallic 
properties, often referred to as "non-Fermi-liquid" (nFL) 
behavior, has triggered a growing interest in quantum- 
critical points (QCPs) that provide a route towards nFL 
behavior [l|. The quantum-critical fluctuations induced 
by suppression of antiferromagnetic (AF) order [2] are 
suggested to mediate the Cooper pairing in exotic super- 
conductors, in an analogous way to phonons in conven- 
tional superconductors [3|. 

Newly discovered iron-pnictide high- Tc superconduc- 
tivity also appears where antiferromagnetism is sup- 
pressed via chemical substitution or pressure [J| . The ex- 
istence of a QCP in iron pnictidcs has been suggested ^- 
|9|. Spin-fluctuation- mediated superconductivity associ- 
ated with the suppression of the antiferromagnetism is 
one likely scenario [l3| , but the identification of the mech- 
anism is far from settled [ll|, [l2| . The difficulty in exam- 
ining the superconducting (SC) mechanism could arise 
from complexity in the materials that can lead to am- 
biguous interpretations; e.g., non-universal SC gap func- 
tions and limited sample quality. It is thus essential to 
find a suitable model system to examine the mechanism 
of superconductivity. 

The isovalent-doped BaFe2(Asi_a;Pj;)2 can be used as 
such a model system. It has the highest Tc (31 K) 
among iron-pnictide superconductors known to have line 
nodes in the SC gap |13l - [la |. Clarifying the mechanism 
that produces its high- Tc nodal gap is thus very impor- 
tant. Since isovalent P-doping is not expected to add 
carriers [1^1, BaFe2(Asi_a;P3;)2 maintains the compensa- 
tion condition, i.e. the volume of the hole Fermi sur- 



faces (FSs) is equal to that of the electron FSs. Very 
clean single-crystals of BaFe2(Asi_a;Pa;)2 allow de Haas- 
van Alphcn (dHvA) experiments that are the most pre- 
cise technique to determine FSs, revealing the detailed 
electronic structure for comparison with band calcula- 
tions [13, [l3l ■ The quasiparticle effective mass increases 
towards the maximum Tc, signaling the enhancement of 
electron-electron correlation. Such an increase in the 
quasiparticle mass as well as nFL behavior inferred from 
resistivity measurements [la, Q^ can be expected when 
the system is in proximity to a QCP. However, direct 
evidence for the existence of a QCP remains lacking in 
BaFe2(Asi_:^P:,)2. 

Here, we report the normal-state spin dynamics in 
BaFe2(Asi_2,Pi:)2 for 0.2 < x < 0.64 investigated by 
■^^P NMR measurements. NMR is highly sensitive to 
low-energy spin fiuctuations and can give information 
of the dynamical magnetic susceptibility. Samples of 
a mosaic of single crystals were prepared as described 
elsewhere [ig]. The dHvA experiments observed sig- 
nals in single crystalline samples from the same batch 
as X =0.41, 0.56, and 0.64, indicating the excellent qual- 
ity of our samples [l7| . The high quality of our samples 
is also suggested by the sharp SC transitions [l6|- '^^P- 
NMR spectra (^^7/27r = 17.237 MHz/T) were obtained 
by sweeping frequency in a fixed magnetic field of 4.12 T. 
The Knight shift K was determined with respect to the 
reference material H3PO4. The '^^P nuclear spin-lattice 
relaxation rate Tj~ was determined by fitting the time 
dependence of spin-echo intensity after saturation of nu- 
clear magnetization to a theoretical 7=1/2 curve with 
a single component of Ti . 

Figure 1 (a) displays '^^P NMR spectra obtained from 
the BaFe2(Asi_a;P2,)2 crystals. Each ^^P NMR spectrum 
consists of a single line, ruling out microscopic inhomo- 
gcneity caused by P substitution. K, which measures 
the effective field at the nucleus produced by electrons. 



is described a.s K = Jv, 



^ chcm 



i^spin is the spin part 
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FIG. 1: (Color online) P substitution evolution of (a) ^^P 
NMR spectra and (b) ^^P Knight shift C^K) determined at 
the spectral peak, obtained in a mosaic of the single crystals at 
4.12 T. (b) The arrows indicate the chemical shift -R'chem (see 
text). The abrupt decrease in K at low temperatures is due 
to the onset of superconductivity. Inset: "^^K vs {TiT)~^'^ at 
270 K for different P concentrations x. 

of K and is related to the uniform spin susceptibility 
x(q = 0), which is proportional to the density of states 
at the Fermi energy N{Ef)- -K'chcm is the chemical shift, 
which is unrelated to x{q = 0)i ^nd is estimated to be 
^ 0.018 % as follows. Since no obvious AF fluctuations 
were detected by NMR at high temperatures as seen in 
Fig. 2, it would be a good approximation to assume that 
{TiT)~-^/^ is proportional to N{Ef) at high tempera- 
tures, i.e. to assume that the usual Korringa relation 
holds at 270 K. Based on the plot of (TiT)-^^^ against 
iiT at 270 K for different x shown in the inset of Fig. [I] we 
can estimate ii'chcm as the intercept. The obtained -ftTchcm 
is 0.018±0.019 %, indicating that i^spin accounts for 86 % 
of the observed Knight shift for x = 0.33. Note that this 
-K^chcm would be a reasonable value, since the chemical 
shift for ^^P in many diamagnetic insulators is of the or- 
der of some hundreds of ppm, which is comparable to this 
Kchcm HI- By assuming K^pin oc x(q = 0) = h%N{Ef), 
the P-substitution dependence of ifspin at 270 K sug- 
gests that the change in N{Ef) would be at most 10 % 
for X < 0.64, which is quantitatively consistent with the 
result of our band calculation discussed below. 

Band-structure calculations by local-density- 
approximation (LDA) were performed for non-spin- 
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FIG. 2: (Color online) The '^^P nuclear spin-lattice relaxation 
rate divided by temperature {TiT)~^ for BaFe2(Asi_i,Pa;)2 
at 4.12 T. Solid lines represent fits to {TiT)~^ = a + b{T + 
9)~^ (see text). As AF fluctuations are suppressed as inferred 
from the suppression of (TiT)~^, Tc (denoted by arrows) also 
decreases. Inset: Fitting parameters of (TiT)"^. The fitting 
parameters a and b are plotted against P concentration x. 
The a and b weakly depend on x, but 9 shows a strong x 
dependence (see Fig. [3|. The small value oi 8 at x = 0.33 
{9 ~ 0) are insensitive to the fitting parameters of a and b. 



polarized BaFe2As2 and BaFe2P2, using the WIEN2k 
package in the APW-f local orbital basis [21[. In addi- 
tion, to obtain systematic changes of the electronic band 
structure for BaFe2(Asi_3;Pa;)2, we performed the LDA 
calculations for three virtual materials with linearly- 
interpolated z = zo{l — x) + zix; (i) BaFc2As2 with fixed 
{a,c) = (ao,Co), (ii) BaFe2As2 with linearly-interpolated 
(a, c) = (ao(l — x) + aix, co(l — x) + cix), (iii) BaFe2P2 
with fixed {a,c) ~ (ai,ci), where ao,i, co,i, and zo,i are 
the experimental values for the crystallographic param- 
eters of BaFc2As2 and BaFc2P2, respectively [la, [23 |. 
N{Ef) barely changes for x < 0.5, and then decreases 
for x > 0.5 (see supplementary information). Such 
behavior is consistent with our Knight shift results. 

As shown in Fig[l] (b), K is almost temperature inde- 
pendent for X < 0.56; and that the absolute value of K 
in the normal state decreases only slightly upon P sub- 
stitution. These data indicate that P substitution does 
not produce significant changes in x(q = 0) and N{Ef)- 
This is in stark contrast to carrier-doped iron-pnictide 
superconductors; in electron-doped Ba(Fei_yCoj,)2As2, 
the Knight shift data indicate that N{Ef) of non-SC 
y — 0.26 is approximately 50 % that of y = 0.08 with 
the maximum Tc of 26 K Q. Such drastic effects on 
N{Ef) via electron-doping is expected from the charac- 
teristic band structure [231; the calculated N{Ef) rapidly 
changes near Ef with a negative gradient, resulting in a 
rapid decrease of N{Ef) with electron-doping. Large 
changes in FSs via charge-carrier doping thus necessarily 
involve dramatic modification in the N(Ef) and Fermi- 
surface nesting resulting in changes in spin fiuctuations. 
In addition to possible changes in Tc due to the modifi- 



cation of spin excitation spectrum \& , drastic changes in 
N{Ep) can also affect severely Tj, [2J|, and the suggested 
giant magnetoelastic coupling may lead to further sup- 
pressions of Tc [23- Therefore, the decrease in N{Ef) 
as well as the suppression of spin fluctuations should 
be taken into account for the interpretation of possible 
changes in Tc for electron-doped Ba(Fei_yCoj,)2As2. In 
contrast, the nearly unperturbed Kgpin by isovalent P- 
doping demonstrates that BaFe2(Asi_a;Pa;)2 is an ideal 
model system to test the relevance of spin fluctuations to 
superconductivity. 

Significant low-energy AF fluctuations arc probed 
near the maximum Tc via Tj~ . {TiT)~^ is described 
by the wave-vector average of the imaginary part of 
the dynamical susceptibility x"{'ii<^o), i-C-, {TiT)~^ ex 
^ \A{q)\'^x"{lj'^o)/^o where ^(q) represents the hy- 
perfine coupling between ^^P nuclear spins and the sur- 
rounding electrons, and wq is NMR frequency. The Kor- 
ringa law TiTK^ = const, generally holds in a FL state, 
but is inapplicable to BaFe2(Asi_a;Pz:)2 near the AF 
phase since Curie- Weiss (CW) behavior is observed in 
(TiT)-! (see Fig. [g. For x = 0.33, (TiT)-^ increases 
significantly down to Tc whereas the Knight shift is con- 
stant. These (TiT)~^ and K data demonstrate convinc- 
ingly that AF fluctuations with flnite q continue to grow 
down to Tc at optimal doping |14| . 

The AF fluctuations in BaFc2(Asi_2;Pa;)2 are enhanced 
signiflcantly as the P concentration is reduced towards 
the maximum Tc (x ~ 0.33), as evidenced by the rapid 
increase in (TiT)^^ from conventional FL behavior at x 
= 0.64 (where (TiT)^^ and K arc almost constant). The 
crossover from FL to CW behavior in {TiT)~^ correlates 
perfectly with the change in the resistivity results [16|; 
As the system evolves from a Fermi liquid {x = 0.71) 
towards the maximum Tc (x = 0.33) near the AF phase, 
the temperature dependence of the resistivity changes 
from T^ to T-linear, one hallmark of nFL behavior. The 
exponent of the temperature dependence of the resistiv- 
ity is shown as a contour plot in Fig. [3] (a). Speciflcally, 
the CW behavior of {TiT)~-^ and the T-linear resistivity 
at x = 0.33 can be explained by the existence of two- 
dimensional (2D) AF spin fluctuations in the theory of 
nearly AF metals [2g| . Indeed, the evolution of the AF 
spin excitations measured by (TiT)~^ upon P substitu- 
tion can be fit consistently by the equation expected from 
the same theory ^, [T^T)-^ = a + b{T + 0)-'^ (see solid 
lines in Fig. [2]). Such 2D AF fluctuations were experi- 
mentally observed in t he p arent BaFe2As2 via neutron 
scatte ring experiments |23l. According to band calcula- 
tions [13J, substantial 2D AF fluctuations can be gen- 



erated by the inter-band spin excitations between the 
multiple FSs predominantly derived from Fe d electrons. 
Our central finding in BaFe2(Asi_2;P2;)2 is that the 
2D AF fluctuations of a quantum-critical nature have 
a clear correlation with the enhancement of quasipar- 
ticle effective mass and Tc as summarized in Fig. [3l 
We found that the Weiss temperature 9 obtained from 
the fltting increases with P substitution and becomes 



almost zero near x ~ 0.33 where the maximum Tc is 
achieved. 9 = K implies that the dynamical suscep- 
tibility probed by {TiT)~^ measurement diverges at ab- 
solute zero, or that the magnetic correlation length con- 
tinues to increase down to T = K. Our result thus 
strongly suggests the presence of an AF QCP near the 
maximum Tc in proximity to the AF phase boundary. 
As the P concentration is varied towards optimal doping 
{x ^ 0.33) from x = 0.64 where the FL state is observed, 
the magnetic fluctuations become dramatically enhanced 
as 9 decreases. Importantly, the quasiparticle mass and 
Tc increase as 9 approaches K. This strongly suggests 
that the AF quantum-critical fluctuations lead to strong 
mass renormalization and unconventional superconduc- 
tivity. Our systematic NMR measurements, which are 
compared with transport measurements, thus provide the 
flrst evidence that the quasiparticle mass enhancement is 
strongly coupled to the AF quantum-critical fluctuations 
in iron-pnictide superconductors as previously observed 
in heavy fermion systems [25 . Furthermore, since N{Ef) 
generally correlates with Tc in conventional BCS super- 
conductors, the enhancement of Tc with approaching the 
QCP from the overdopcd side cannot be accounted for 
by an nearly unperturbed N(Ep), demonstrating clearly 
that superconductivity in BaFe2(Asi_a;Pa;)2 is tuned pre- 
dominantly by the AF fluctuations. Therefore, we con- 
clude that the 2D AF quantum-critical fluctuations are 
likely to play a central role in the occurrence of uncon- 
ventional superconductivity in BaFe2(Asi_2;P2-)2- 

Previous NMR measurements also indicated that there 
is a strong connection between low-energy spin fluctu- 
ations and superconductivity in 122 [ij]. 111 [20| and 
11 iron-based superconductors [301. These suggest that 
the relevance of spin fluctuations to the high- Tc super- 
conductivity is universal despite the non-universal gap 
structure. The strong connection would be naturally un- 
derstood if one considers the superconductivity observed 
near magnetic phase and the well-nested hole and elec- 
tron Fermi surfaces. By contrast, we reported that the 
connection is relatively weaker in LaFeAs(Oi_a;Fa;) [3l| . 
Such a weak link may originate from much lower Tc of 
26 K in LaFeAs(Oi_:rF:r) than that of i?FeAs(Oi_;j;F^) 
(R = Ce, Pr, Sm, etc.). However, it remains unresolved 
whether spin fluctuations are also important for high-Tc 
superconductivity in the 1111 system with Tc exceeding 
50 K, since magnetic rare-earth atoms R hinder quasipar- 
ticle excitations responsible for superconductivity j32i |. 
Therefore, examinations of the connection between mag- 
netic fluctuations and superconductivity of high-Tc 1111 
iron pnictides would be crucial for establishing general 
view of iron pnictide superconductivity. 
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FIG. 3: (Color online) Temperature - P concentration phase 
diagram of BaFe2(Asi_a;Pa:)2. Tc, Ts, and Tn denote the 
SC, orthorhombic-tetragonal, and AF transition tempera- 
tures, respectively, (a) Colors represent the exponent, a = 
d(ln Apa6)/d(lnr), at zero magnetic field for the resistivity 
curves in Ref. 24, where Apab = pab — pab{T — )■ K). (b) The 
triangles represent 9 obtained by fitting (TiT)^^ data to the 
equation, {TiT)~^ = a-\-h{T + 6)~^ (error bars are within the 
symbol size). T- linear resistivity, suggesting quantum-critical 
behavior, is observed in proximity to a magnetic instability 
signaled by S ~ at a; = 0.33. These quantum-critical fluctu- 
ations are found to correlate with the enhancement of the ef- 
fective mass as denoted by the pentagons (taken from Ref. 17). 
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